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We study the diffusion of tracers (self-diffusion) in a homogeneously cooling gas of dissipative 
particles, using the Green-Kubo relation and the Chapman-Enskog approach. The dissipative par- 
ticle collisions are described by the coefficient of restitution e which for realistic material properties 
depends on the impact velocity. First, we consider self-diffusion using a constant coefficient of resti- 
tution, e =const, as frequently used to simplify the analysis. Second, self-diffusion is studied for a 
simplified (stepwise) dependence of e on the impact velocity. Finally, diffusion is considered for gases 
of realistic viscoelastic particles. We find that for e =const both methods lead to the same result 
for the self-diffusion coefficient. For the case of impact-velocity dependent coefficients of restitution, 
the Green-Kubo method is, however, either restrictive or too complicated for practical application, 
therefore we compute the diffusion coefficient using the Chapman-Enskog method. We conclude that 
in application to granular gases, the Chapman-Enskog approach is preferable for deriving kinetic 
coefficients. 



The diffusion of tracer particles (self-diffusion) 
in a force-free gas of dissipatively colliding par- 
ticles, also called granular gases, is studied. 
Main theoretical approaches for the calculation 
of transport coefficients are the Green-Kubo and 
Chapman-Enskog methods. The former is based 
on fluctuation-dissipation relations, which ex- 
press transport coefficients in terms of the time 
correlation function of dynamical variables. The 
latter approach is based on the Boltzmann equa- 
tion, that describes the evolution of the distribu- 
tion function of tracers. The dissipative nature 
of the gas is characterized by the coefficient of 
restitution, e which is the central quantity in the 
theory of granular gases. This coefficient quanti- 
fies the loss of kinetic energy due to a single par- 
ticle collision. We consider three different mod- 
els for the coefficient of restitution: (i) e=const., 
(ii) a stepwise dependence of e on the impact ve- 
locity, which mimics the basic property of this 
coefficient, and (iii) a realistic model for e as 
a function of the impact velocity, derived from 
the contact mechanics of viscoelastic bodies. For 
the simplest model (i), both theoretical methods 
yield the same coefficient of diffusion. For the 
more realistic models (ii) and (iii), however, the 
Green-Kubo approach is either restrictive or too 
complicated for practical application. Therefore 
we conclude that the Chapman-Enskog method is 
preferable for deriving kinetic coefficients of gran- 
ular gases. 



I. INTRODUCTION 

A homogeneously initialized force free granular gas, 
i.e., a rarefied systems of inelastically colliding macro- 
scopic particles, stays homogeneous during the first stage 
of its evolution while its temperature, that is the aver- 
age kinetic energy of particles, decays due to dissipative 
collisions. This state of a granular gas is called the ho- 
mogeneous cooling state. At later stages, granular gases 
develop density inhomogeneities 0, and pronounced 
spatial correlations of the velocity field Throughout 
this article we consider granular gases in the homoge- 
neous cooling state. 

Similar as ordinary molecular gases, granular gases 
in the homogeneous cooling state are characterized by 
their temperature and their velocity distribution func- 
tion, which is close to the Maxwell distribution. Like in 
molecular gases, impurities (guest particles), that is par- 
ticles which differ from the gas particles by mass, size or 
other parameters, move irregularly due to random colli- 
sions with the surrounding gas particles. This Brownian 
motion of guest particles, i.e. diffusion, leads to uni- 
form spreading of an initially localized ensemble of im- 
purities over the available volume, ff the guest particles 
are mechanically identical to the host particles, but may 
be somehow distinguished (e.g. by color), the process 
is called self-diffusion and the guest particles are called 
tracers. 

The processes of diffusion and self-diffusion in granu- 
lar gases differ from those in molecular gases because of 
the dissipative nature of particle collisions. A collisions 
of dissipatively colliding particles is characterized by the 
coefficient of restitution e. This coefficient relates the 
post-collision velocities v- and v- of the colliding parti- 
cles i and j to their pre-collision velocities u, and Vj . For 
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particles of identical mass the collision law reads 
v t =v l — {Vij ■ e ) e 

Vj = Vj + i v ij ■ e ) e » 



(1) 



where Vij = Hi — 
fij = fi — fj describes the relative particle position at 
the instant of the collision. For particles of different mass 
the corresponding relation has a similar form jj|. The 
collision law Q takes only the normal component of the 
particles' relative velocity into account. Correspondingly, 
e, is also called the normal coefficient of restitution - 
the tangential component of the relative velocity is not 
affected by a collision (see Q, Chap. 3.4 for a detailed 
discussion of this non-trivial approximation). In what 
follows wc consider the normal coefficient of restitution. 

Granular gases are never in thermodynamic equi- 
librium because of the permanent loss of energy due 
to dissipative particle collisions. Therefore, the basic 
law of equilibrium statistical mechanics, the cquipar- 
tition of kinetic energy, is not applicable for granular 
gases. Consequently, the temperature of the gas parti- 
cles may differ from the temperature of t he g uest par- 
ticles j||lllllll|l3|ll|ll|ll|ll|li. Hence, 
in granular gases the temperature is species dependent. 
The violation of equipartition leads to many interest- 
ing phenomena, like criticality of the Brownian motion 
|14L Il5| and turns the analysis technically complicated. 
Here we restrict ourselves to the case of self-diffusion, 
when the gas may be characterized by a single tempera- 
ture, that is, we discuss the diffusion of tracers. In spite 
of being much simpler, self-diffusion demonstrates the 
most salient properties of transport in granular gases. 
It allows to illustrate the application of the main theo- 
retical approaches for the computation of the transport 
coefficients, the Green-Kubo and Chapman-Enskog ap- 
proaches, without much technical complications. 

We show that for the case of a constant coefficient of 
restitution, the Green-Kubo approach, which has been 
initially elaborated for equilibrium systems, may be ef- 
ficiently exploited. For a direct and thus simple appli- 
cation of this method, the dissipative gas dynamics can 
be mapped onto the corresponding dynamics of a system 
which remains in an effective steady-state. We will show 
that in this case the Green-Kubo and Chapman-Enskog 
methods yield the same result for the self-diffusion coef- 
ficient. For the case e =const. such a simple map exists. 

The assumption of a constant coefficient of restitution 
simplifies the analysis drastically, however, it does not 
agree with experimental observations fl7L llSj . More- 
over, the assumption e =const. contradicts the mechan- 
ics of materials [l9l |2(1 12 ll \'2'A l23l l24j and even violates a 
dimension analysis [25j. Instead, experiments as well as 
theory show that the coefficient of restitution increases 
with decreasing impact velocity. Thus, for g — > par- 
ticles tend to collide elastically, that is e — > 1 if surface 
effects are neglected. 



This basic property motivates a drastically simplified 
collision model |2al27|: 



<9) 



for g > g* 
for g < g* 



(2) 



The unit vector e = fy/Vy with where g 



e is the normal component of the rela- 



tive velocity at the impact. Thus, it is assumed that 
impacts at large normal relative velocity occur with a 
constant coefficient of restitution £*, while collisions at 
small g are elastic, i.e., the threshold g* separates elas- 
tic from inelastic collisions. The model parameters, e* 
and g* , are neither related to the material properties of 
particles nor derived from first principles of the collision 
dynamics. 

The coefficient of restitution as a function of the im- 
pact velocity may be also derived analytically by inte- 
grating Newton's equation of motion for a colliding pair, 
provided the particle interaction force is known. The sim- 
plest physically realistic mechanical model of dissipative 
collisions assumes that particles interact by viscoelastic 
forces 0, 0, HH • If inelasticity is small, the solution 
of the contact problem for viscoelastic spheres yields the 
coefficient of restitution HH Hal 
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(3) 



where the coefficient 7 is a known function of the elastic 
and dissipative constants, i.e. of the Young modulus and 
the Poisson ratio of the particles material, and of the 
masses and radii of the particles (see j^,!^ for details). 

We will show that for the case of impact-velocity de- 
pendent s there is no simple map of the dynamics of the 
dissipative system onto a corresponding dynamics of a 
steady-state system. This makes the application of the 
Green-Kubo method problematic. Therefore, for the case 
of an impact- velocity dependent coefficient of restitution 
the Chapman-Enskog approach is used. We apply the 
latter method to calculate the self-diffusion coefficient 
for the case of a simplified impact-velocity dependent 
coefficient of restitution as well as for realistic gases of 
viscoelastic particles. 

The rest of the paper is organized as follows: In Sec. II 
we discuss the microscopic interpretation of the macro- 
scopic diffusion coefficient which suggests a route of its 
evaluation. The generalization of the diffusion coefficient 
to non-equilibrium systems, such as granular gases is also 
discussed. In Sec. Ill, the formalism of the pseudo- 
Liouville operator is addressed, as used for hard-sphere 
systems, and Sec. IV illustrates its application to gran- 
ular gases in the homogeneous cooling state. Three dif- 
ferent collision models introduced above are there ad- 
dressed. The Green-Kubo approach is discussed in Sees. 
V and VI, while the Chapman-Enskog method is de- 
scribed in Sees. VII and VIII. In Sees. V, VI and VII we 
explain the general concept of these methods for the sim- 
ple case e = const. Then we consider the more general 
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case of impact-velocity dependent coefficient of restitu- 
tion. We summarize our findings in Sec. IX and in Ap- 
pendixes A and B some technical detail are given. 



II. MACROSCOPIC AND MICROSCOPIC 
MEANING OF THE DIFFUSION COEFFICIENT 

Consider a uniform gas of hard spheres of mass m and 
diameter a whose particles collide due to the collision law 
Eq. (JJJ, specified by the coefficient of restitution e. Let 
tracers particles be embedded in the gas and let f(f, v, t) 
and f s (v a ,f s ,t) be respectively the one-particle distribu- 
tion function of the gas particles and of the tracers. Then 
the local densities are defined by 



n s (r s ,t) = Jdv s fM,r s ,t) 
n(f,t) = / dvf(v,f,t) . 



(4) 



If the concentration of the tracers is not uniform, a dif- 
fusion flux of these particles arises, which is directed op- 
posite to the concentration gradient: 



J s (r ) = -DVn s (r ; 
With the continuity equation 

dn s (f ) 



dt 



V J s (f ) = 



we obtain the diffusion equation 



dn s (f) 
dt 



DV 2 n s <r' 



(5) 



(6) 



(7) 



The diffusion coefficient D in Eq. J£J, thus, relates the 
macroscopic flux of particles J s (r ) to the macroscopic 
concentration gradient Vrt s (f ). 

However, there is also a microscopic interpretation of 
the diffusion coefficient: It defines the time dependence 
of the mean square displacement of a tracer particle. As- 
sume at t — there is an ensemble of impurities located 
at the origin, r = 0. Then the mean square displacement 
of the tracers at time t reads 



([r s (t)] 2 )^ f df[f(t)] 2 n s (r,t) 



(8) 



where N s — J dfn s (r, t) = const, is the total number of 
tracers. Now we multiply both sides of Eq. J7J by f 2 /N s 
and integrate over df : 



d_J_ 

dtN~ s 



r— / dfr 2 n s = D j drr 2 V 2 n s . 



(9) 



which turns by means of Eq. JSJ and Green's theorem 
|2Sj for the right-hand side into 

j t ([r s (t)} 2 ) =D±[ drn s V 2 r 2 = 6D , (10) 



that is, we obtain a relation between the mean square dis- 
placement of the tracers and the coefficient of diffusion: 



[r s (t)] 2 ) =6Dt. 



(11) 



This equation relates the diffusion coefficient D, i.e. a 
macroscopic quantity, to a microscopic quantity, namely 
the mean square displacement. Therefore, Eq. (|II|) sug- 
gests a method to compute the diffusion coefficient by 
evaluating ([f s {t)] 2 ^. To find this quantity we apply the 
kinematic relation 



f s (t) = / v s (ii) dtx (f s = at t = 0) (12) 
Jo 



and write 



([r s (i)] 2 ) = Q v s {ti)dh ■ j v s (t 2 )dt 2 \ . (13) 

Using the symmetry of the velocity autocorrelation func- 
tion, 

K v (ii,i 2 ) = (v s (tr) • v s (i a )) = K v (t 2 ,h) , (14) 
we can recast Eq. p3|l into the form 

([f s (0] 2 ) = 2^ dhj^ dtaK v (t a ,ti) . (15) 

If the system is in a steady state its properties are invari- 
ant under a shift of time, thus, the velocity autocorrela- 
tion function depends only on the time lag r = t 2 — t\, 
that is, K v (t 2 ,ti) — K v (t) and decays with a charac- 
teristic relaxation time t v . Transforming the variables 
(ti, t 2 ) — > (r, ti) with n = t\ the mean square displace- 
ment reads 



([r s (if) = 2t j Q dr (v s (0) • v s (r)) (l-j). (16) 

For large time as compared with the characteristic time, 
i.e. for t 3> t v , we obtain 

([r 3 (<)] 2 ) =2t j dr (v s (0) • v s (r)) = 6Dt (17) 



and consequently the diffusion coefficient 
D =ll (v s (0)-v s (t))dt. 



(18) 



Equation HI8|I expresses the macroscopic transport co- 
efficient D by the time integral of the microscopic ve- 
locity correlation function. This fluctuation-dissipation 
relation expresses a dissipative quantity (diffusion coeffi- 
cient) via a fluctuating quantity (particle velocity). Sim- 
ilar expressions relate other transport coefficients such as 
viscosity, thermal conductivity, etc. to the corresponding 
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fluctuating quantities [29j. The method of calculation 
of the transport coefficients via the microscopic time- 
correlation functions is called Green-Kubo approach. Ini- 
tially, the Green-Kubo theory has been developed for 
equilibrium systems of elastic particles (e.g. |29fl). Later 
it has been generalized for dissipative gases [20> El • 
For the case of self-diffusion the velocity correlation func- 
tion in Eqs. (|13I18|) for the tracer particles is equal to 
that of the surrounding gas particles. 

The derivation given above applies to equilibrium sys- 
tems or systems in a steady state. Granular gases are, 
however, never in a steady state due to persistent cool- 
ing. Therefore, the concept of the diffusion coefficient 
has to be generalized to non-equilibrium systems and a 
time-dependent diffusion coefficient D{t) (also called dif- 
fusivity) must be used |33ll34j : 



{[f{t)f) = I D{t')dt>. 



(19) 



As we show below, in certain cases, the dynamics of gran- 
ular gases may be mapped onto the dynamics of a steady- 
state system. In this case the above Green-Kubo deriva- 
tion of the diffusion coefficient is still applicable. 



Vij = iJi — Vj and e = fj/fij. The factor \vij ■ e | in the 
integrand of Eq. I|22l) gives the length of the collision 



cylinder, the Hcavisidc step function (- 



e ) selects 



approaching particles and the S function specifies the unit 
vector e at the instant of the collision. The operator 
is defined by 



b ij F (ri,rj,Vi,Vj 



) = F(? i ,r j ,vl,v' j ---) , (23) 



where F is some function of the positions and coordi- 
nates. The post-collision velocities of the colliding pair, 
v[ and v(, are given in terms of their pre-collision val- 
ues according to Eq. . The pseudo-Liouville operator 
defined above describes the evolution of dynamical vari- 
ables, such as positions of particles, their velocities, etc. 
It has been introduced first to describe the evolution of 
a system of elastic hard spheres [3(| |37j ■ 

The velocity time correlation function may be found 
by performing the formal integration of Eq. (|20|l for the 
variable A = v for t > t': 



v(t) = e c ( t - t ')v(t>) 



(24) 



yielding 



III. DYNAMICS OF DISSIPATIVE HARD 
SPHERE SYSTEMS: GENERAL APPROACH 

Because of the singular character of the hard-sphere 
interaction (the interaction potential jumps to infinity at 
the point of particles contact) it is not possible to apply 
directly Hamilton's equations to describe the dynamics of 
hard-sphere systems. Nevertheless, the dynamics is still 
completely deterministic. Indeed, between collisions the 
particles move along straight lines at constant velocities, 
and the post-collision velocities are uniquely determined 
by the collision law Eq. Q). This allows to construct 
the pseudo-Liouville operator C which governs the dy- 
namics of the system. Let A(t) be an arbitrary function 
of particle coordinates and velocities, then the evolution 
equation for this quantity reads 



d 
dt 



A(fv) = CA(r,v) , 



(20) 



where the pseudo-Liouville operator has the form (e.g. 



N N N 

i— 1 j>i 



(21) 



i=l 



with the binary collision operator 



deO {—Vij ■ e) \fij ■ e \ 6 (fj — ae ) 



Here N is the number of particles, fi and Vi are respec- 
tively the position and velocity of particle i, fij = fi — fj, 



(v(t')v{t))= f dTp{t')v{t')e c ( t - t ')v(t') 



(25) 



where J dT denotes integration over all degrees of free- 
dom and p(t) = p(fi, . . . ,fisr,vi, . . . ,VN,t) is the N- 
particle distribution function which depends on temper- 
ature T, particle number density n, etc. Although Eq. 
(|25|l provides the velocity correlation function K v , whose 
time integral yields the diffusion coefficient, it may not 
be easily used in this general form. Instead, for prac- 
tical application of the Green-Kubo method, in general, 
further approximations are required. 



IV. HOMOGENEOUS COOLING STATE 
A. General properties 

We consider a granular gas in the homogeneous cooling 
state when the number density of the gas n = N/V is 
uniform in the gas volume V, while the gas cools down 
due to inelastic particle collisions. 

The evolution of the gas temperature may be easily 
obtained from its definition as the mean kinetic energy 
of the particles, 



—nT(t) — n ( —mv 2 
2 w \2 



-f(f,t)dv. (26) 



From the definition of the Liouville operator, Eqs. I|21t - 
1221 follows 



3 dT 
2 ~dt 



m d i 9 \ 
25^ 



rn 
T 



(27) 
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where (. . .) t denotes averaging using the iV-particle dis- 
tribution function p(t) at time t. With Eqs. H21I22[1 and 
exploiting the identity of the particles we obtain 



3 dT 
2 H 



J2 



(N 



(28) 



Using the definition of the binary collision operator, Eq. 
, and 



(8& - 1) «? = K) a - 

Eq. (|2H1) turns into 
3 dT 



(29) 



2 eft 



(N — 1) / dri...dr N j dv% . . . di)Np(t)x 



deQ (—vi2- e)\vi2- e\ 5 (ri2 — ae) (v[) 



(30) 



Usin g th e definition of the two-particle correlation func- 
tion |23 and the hypothesis of molecular chaos, for a 
uniform system we may further write 



N(N - 1) J df 3 ... df N J dv 3 ... dv N p{t) 

= h {ri,r 2 ,v u v 2 ,t) ( 31 ) 

= 92 (ru) f (vi,t) f {v 2 ,t) , 

where g 2 (^12) is the pair correlation function. Substitut- 
ing finally the latter expression into Eq. (fTTTTfl yields 



3 dT 
2 ~dt 



m 
2ri 



92 (ct) <j 2 J dvidv 2 J de<d (-vi 2 ■ e] 



(32) 



Introducing the time-dependent thermal velocity, the 
reduced velocities, and the scaled velocity distribution 
function defined respectively by 



v T {t) 



f(v,t) 



4/(c,t) (33) 



we write the last equation in a more traditional form (e.g. 



1) 



C = ^na 2 g 2 (a) v T p 2 



(34) 
(35) 

p 2 = -i J dc\dc 2 JdeQ(-ci 2 -e) (36) 
x |c 12 ■ e I / t) f (c 2 , t) A (c? + c|) , 

where the symmetry with respect to the exchange of par- 
ticle indices 1 «-> 2 has been exploited. The notation 



A (ip) = ip' - ip 



(37) 



describes the change of a dynamical variable ip due to 
a collision, with ip' being the after-collision value. The 
contact value of the pair correlation function reads 



92 (ct) 



(1-77)3' 



1 3 

with 7] = —nna . 
6 



(38) 



For a granular gas in the homogeneous cooling state 
the reduced distribution function / (c, t) is close to the 
Maxwellian 0, IH IH E E3 ■ It may be written with 
good accuracy, keeping only the first two non- vanis hing 
terms of the Sonine polynomial expansion 0, 13^. l4pl l4l| : 



/(<?,*) 



[l + a 2 (t)S 2 (c 2 )] 



where 



(c) =7T- 3 / 2 CXp(-C 2 ) 



is the rescaled Maxwell distribution and 



S2(C 2 ) 



(39) 



(40) 



(41) 



is the second-order Sonine polynomial. The second So- 
nine coefficient a 2 characterizes the shape of the veloc- 
ity distribution function. It sensitively depends on the 
particular collision model, which we discuss separately 
below. 

According to Eq. (|35|l . to find the cooling coefficient 
£ we need to compute p 2l which may be written as an 
integral of the general form 

jdcx Jdc 2 JdeO {-c 12 ■ e) \c 12 ■ e | </>(ci)0(c 2 )Expr (■ • ■ ) , 

(42) 

where 



Expr= --A(c 2 



[1 + a 2 S 2 (c 2 )] [I + a 2 S 2 (c 2 2 )] 



(43) 

We call integrals of the type Eq. I|42|) kinetic integrals 
0, |4l|. The kernel Expr (...) may contain the pre- and 
post-collision velocities, c\/ 2 , Cy 2 , the Sonine polynomi- 
als Si of these velocities, the Sonine coefficients aj, the 
unit vector e, the relative velocity c*i 2 , the center of mass 
velocity of the colliding particles, C, scalar products of 
the mentioned vectors, and other variables. Further, the 
coefficient of restitution may enter the kernel Expr (...), 
either e = const or as a function of the impact velocity 
e — e(g) as discussed below. 

The methods for co mput ing kinetic integrals are ex- 
plained in detail in 0, |4l|. Kinetic integrals may be 
represented by sums of standard-type integrals which al- 
low for a solution by means of computational symbolic 
algebra. Many physical quantities of interest may be ex- 
pressed in terms of kinetic integrals 0, . Let us con- 
sider in more detail the collision models discussed in the 
Introduction. 
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B. Granular gases with e=const 

As mentioned above, the coefficient of restitution is 
a function of the impact velocity. Nevertheless, here we 
assume e =const. to simplify the explanation of the basic 
approaches. For this case the second Sonine coefficient 
reads 0, 0] ( see Appendix [XJ where the derivation of 
0,2 is outlined): 



16(l-e)(l-2e 2 ) 
81- 17e + 30e 2 (l-e) 



(44) 



Equation l|44J) gives the second Sonine coefficient in a 
linear approximation with respect to 02, the next-order 
solution is derived in [l^. Note that for e=const. the 
shape of the scaled distribution function / and, hence, 
the Sonine coefficient do not depend on time. 1 

The coefficient /12 may be found by means of the com- 
putation method elaborated for the kinetic integrals 01 : 



M2 



^(l 



32 



n 2 



(45) 



Substituting /i2 into Eq. i|35|) yields the cooling coeffi- 
cient (,. Since /12 does not depend on temperature for 
e = const., the cooling coefficient scales as £ cx T 1 / 2 , see 
Eqs. (|33I35|1 . Hence, from Eqs. JM] |23J) we obtain for 
the evolution of temperature 



T = 



T(0) 



/'2 



(1+i/ro) 2 



=Tc(0)- 



(46) 



where /i 2 and r c are given respectively by Eq. (|45|) and 
Eq. (53) . 



C. Granular gases with stepwise impact- velocity 
dependent coefficient of restitution 

For granular gases of particles which collide with the 
simplified stepwise impact- velocity dependent coefficient 
due to Eq. (0, the velocity distribution function is 
close to a Maxwell distribution with good accuracy [26| . 
Therefore, a 2 ~ for these gases, and the coefficient /X2 
reads (see Appendix lB|l 



2tt 1 



1 



mg 



4T 



exp 



mg 



4T 



(47) 

Substituting (|47|l into Eq. (|35|l yields the cooling coeffi- 
cient C for this collision model. According to Eqs. 147|) 



1 MD simulations of two-dimensional granular gases show that 
the second Sonine coefficient ai re laxes to the theoretical value 
within a few collisions per particle [43, 44] . At later times, how- 
ever, 02 deviates from the theoretical result due to accumulating 
correlations between the particles |44| . see Appendix 1X1 



and l|35[) , £ decays exponentially fast as a function of tem- 
perature, £ oc exp(— const. /T), that is, the cooling of the 
gas slows down exponentially as the gas evolves. Asymp- 
totically, for t — > 00, this leads to a logarithmically slow 
decay of the temperature j2||: 



T(t) = 9 f , a = 2g (l - £ * 2 ) g 2 (a)a 2 nJ ^ , 
2 log erf V 2m 



where go = 3*/^r(0). 



(48) 



D. Granular gases of viscoelastic particles 

For gases of viscoelastic particles the shape of the ve- 
locity distribution function and, thus, the Sonine coeffi- 
cients depend on time. The time scale of diffusion pro- 
cesses - the hydrodynamic time scale - is much larger 
than the mean collision time. In this case one obtains 
approximate expressions |45j : 



ft 2 = ft2i<5' + a 2 2^ 



12 



(49) 



where 021 ~ — 0.388 2 and 022 ~ 2.752 are pure numbers 
and 5' is the time-dependent dissipation coefficient, 



S'(t) 



J_ 



v T {t)} 



1/5 



(50) 



with Ci « 1.1534. 

Similarly, solving the kinetic integral, the coefficient \i 2 
is obtained for gases of viscoelastic particles, 



^2 = Uq5' - lu 2 5 



12 



(51) 



where u>q « 6.485 and u>2 ~ 9. 

By means of Eqs. I|35|) . (|50|l . and (|51|l we obtain in 
leading-order, ( ~ T 3 / 5 and find 0, ^ 



T 



T(0) 



16 7 /T 



(1+t/ro) 5 / 3 



1/10 



-<y„ — ! — I T c(0) 1 
(52) 

where go — 0.173. Hence different collision models pre- 
dict different evolution of temperature, which asymp- 
totics ranges from the power law, T — > t~ a with a = 2 
for e = const, and a = 5/3 for viscoelastic particles, 
to the logarithmic decay, T — ► (logi) -1 for the stepwise 
impact-velocity dependent e. 



V. VELOCITY CORRELATION FUNCTION 
FOR e=const 

Again we consider a granular gas in the homogeneous 
cooling state and assume e=const. Let time be measured 



2 The numerical values of the constan ts 021, 022, Cx, <^0i <^1, and 
(jo are known analytically I4I H I4TI . 
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in units of the mean collision time T c (t), 

dt f* dt 1 



nth 3 



dr = 



Tc(t) ' 



"(*) 



Top) 



(53) 



which follow from the definitions of vt and T, and intro- 
duce the laboratory-time dependent velocity relaxation 
time 



T„(t) = f v T c (t) 



(59) 



- 1 {t) = A^g 2 {a)a 2 n 



T{t) 



(54) 



The dynamics of the system is stationary when replacing 
the true (laboratory) time by the rescaled time r. Ne- 
glecting particle-particle correlations (see Q| ) the dynam- 
ics of the system equals the dynamics of an equilibrium 
system. Thus, by rescaling time, the dynamics of the dis- 
sipative gas is mapped onto the stationary dynamics of 
a gas of elastically colliding particles jl, l^ll E3 • In this 
state the collision frequency and temperature are con- 
stant and the only difference from an equilibrium molec- 
ular gas is that the collisions still follow the rules of the 
dissipative impact, Eq. QJ, with the coefficient of resti- 
tution e. The velocity time-correlation function may be 
then written in the form 

(fli (f) ■ n (t)) = v T (f) v T (t) (c! (O • ci (t)> , (55) 

i.e., we need to find the time-correlation function for the 
reduced velocity c(t). In the new time scale the collision 
frequency and temperature are constant and the system 
is in a stationary steady state. Consequently, the time 
correlation function (ci (r') • c\(t)) docs not depend on 
the times r' and r separately, but only on their difference 
\r'-r\. 

Adopting the molecular chaos hypothesis, the veloci- 
ties of two colliding particles are not correlated. Hence, 
the collision process for a particle is stationary and 
Markovian. If wc additionally assume that the process is 
Gaussian, we deduce from Doob's theorem |2^| an expo- 
nential velocity correlation function: 



( c l( r ')-cl(r))=( C ?> e xp(-L^) 



= -exp 



(56) 



with t v being the velocity relaxation time, measured in 
units of the mean collision time r c (t). Obviously, f v is 
constant since the system is in a stationary state. 

To compute the time derivative of the correlation func- 
tion at time t' = t, we use the relations 



dvT 
~dT 

f«r(*) 



= <«?(*)>, 



-\c(t)v T (t) 



(57) 
(58) 



3 To be precise, r c differs slightly from the actual mean colli- 
sion time due to deviations of the velocity distribution from the 



which corresponds to the reduced relaxation time t v . 
Eventually, we obtain 



dt 



(*(*0 •*(*)>!*=* 



<«?(«)> 



(60) 

On the other hand, using Eq. (|25J) . we obtain for the 
same quantity 



(n(t')-Mt))\ t >=t 

= | f dTp{t')v 1 -e^- t ')v 1 



= (N-l)(vi-f l2 vi} t 



(V! ■ Cvi) 



(61) 



The last equality in Eq. (|61(l follows from the properties 
of the Liouville operator, Eq. (|21l I22II and from the iden- 
tity of the particles. Comparing Eqs. 161(1 and (|60|l we 
obtain the velocity relaxation time r v {t): 



r-\t) = -{N-l) 



(vi ■ T X2 Vi / 



j i it 



(62) 



Maxwell distribution Q. 



Physically, the first term on the right-hand side describes 
the relaxation (i.e. the decay) of the velocity correla- 
tion function due to collisions which randomize the di- 
rections and amplitudes of the particle velocities. It has 
the same nature as for molecular gases where particles 
collide elastically. The second terms describes the decay 
of the velocity correlation function due to the decrease 
of the thermal velocity of the cooling gas. For almost 
elastic particle collisions, e < 1, each particle suffers a 
large number of collisions before the temperature decays 
noticeably. Thus, the velocity relaxation time is domi- 
nated by the first term and the decrease of the thermal 
velocity on the time scale of the velocity relaxation t v 
may be neglected, yielding the adiabatic approximation 
of the velocity relaxation time, t v f=a T v>a d- 

Hence, Eqs. (jSSJ, 123), an d lEU define the veloc- 
ity correlation function of granular gases when e =const 
is assumed. Note that in spite of the molecular chaos as- 
sumption, in the laboratory time the velocity correlation 
function does not decay exponentially as it does for gases 
of elastic particles. 
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VI. GREEN-KUBO DIFFUSION COEFFICIENT 
FOR e=const. 



where we take into account that f v = t v /t c and 
2T/m. Hence, D = T v T/m and, therefore, 



To find the diffusion coefficient it is convenient to use 
the reduced time, measured in collision units and the 
corresponding reduced length 



(67) 



r = r r, r = v T (t) r c (t) . 



(63) 



Then the reduced mean square displacement may be ex- 
pressed as the time integral of the reduced velocity cor- 
relation function (see Eq. j!5|l): 



(Af(r)) 2 ) =2£dn£dr 2 ^exp^-^^ 



(64) 

Consider now the time derivative of the reduced mean 
square displacement at time r>r„, 



To find t v _ a d we evaluate (tii ■ T12V1J using the relation 



(vi ■ T12 vx^j = - (vu ■ T 12 vi^ , (68) 



which follows from the symmetry properties of the binary 
collision operator, Eq. J22J). We also use 



dr 



(Af 2 ) = 3 / dnexp 



3f„ 



1 — exp 



T-Ti 



3t v , 



(65) 



V12 ■ 



y>ij - l) «1 = V12 ■ {v[ - Vi) 



1 2 

= -5 i 1 + £ ) ("12 ' e) , 



(69) 



which reads in unsealed variables 



£<^W>=££<^(r)> = ^. 
_ v 2 T {t)rl{t) o _ T{t) 



3t v = 6^t v = 6D(t) , 



(66) 



which in its turn follows from the definition of the oper- 
ator bfj, Eq. (J2SJ. Using Eq. (jHEJ and substituting Eq. 
B69JI into the definition of T v ,ad in Eq. (|62(l . we obtain 
the velocity relaxation time in adiabatic approximation: 



^IdW = - [<«?>,] 1 (N-l)\ (vu ■ T 12 



3T(t) 



(iV — 1) / dfi . . . dfjv / ■ ■ • dvNP (t) a / de 6 (— «i 2 • e ) |«i 2 • e | 5 (fi 2 — ere ) (1 + e) (v\2 ■ e ) z 



(70) 



This equation may be evaluated in the same way as the cooling coefficient £ in Eqs. I|3UI34|I . The result reads 



Tv,\a (*) = gWT(t).92(ff)fT 2 n / dcidc 2 / de6 (-ci2 • e ) |c 12 • e | / (ci, t) / (c 2 , t) (1 + e) (c 12 • e f , ( 71 ) 



which again is a kinetic integral of the form given in Eq. 
(021. Using Eq. (JSSJ for the distribution function / (c, t), 
we can evaluate this kinetic integral by the scheme ex- 
plained in 0, |4l| and find 



diffusion coefficient: 



D(t) 



4D (t) 



(l + e) 2 (l + ^a 2 )' 



(73) 



where D (t) is the Enskog self-diffusion coefficient, 



which together with Eqs. (|34|l and (|45fl for the cool- 
ing coefficient £ leads to the final expression for the self- 



T(t) 



(74) 
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VII. CHAPMAN-ENSKOG APPROACH FOR 
THE DIFFUSION COEFFICIENT FOR e=const. 

We start from the Boltzmann equation for the distri- 
bution function of tracers f s (v, r, t) |j, |48| 



f s (vi,t) 

= g 2 (<j) a 2 / dv 2 / de 9 (— #12 • e) \v 12 ■ e \ x 



(75) 



X I -f s (V?, t) f «, t) - fs (Vl,t) f (V 2 , t) 

g 2 (<r) /(/,/.), 



where /(/, / s ) abbreviates the collision integral, f(v, t) is 
the velocity distribution function of the gas particles, and 
the index s refers to the tracer particles. In Eq. I|75l) . v[' 
and v 2 denote the velocities of particles after an inverse 
collision, 



v 1 =v 1 — {v 12 ■ e)e 

2e 

v 2 =v 2 + (U12 • e ) e . 



(76) 



The corresponding direct collision according to Eq. QJ, 
thus, turns the initial velocities {v" ,v 2 ) m to (vi,v 2 ). 

The factors in the integrand at the right-hand side 
of Eq. (|75|) have their conventional meaning: \vi 2 ■ e\ 
is the length of the collision cylinder, 9 (—V12 • e ) se- 
lects only approaching particles (since particles depart- 
ing from each other do not collide) and g 2 (a) accounts 
for the increased collision frequency due to the finite di- 
ameter of the particles and the corresponding excluded 
volume. Since the concentration of tracers is small, one 
can assume that they do not affect the velocity distribu- 
tion function of the gas particles, which is given by the 
solution, Eq. for the homogeneous cooling state. 

Solving the Boltzmann equation for f s (r , v, t) , which 
with a given function / (v, t) is also called Boltzmann- 
Lorentz equation, one finds the diffusion flux 



Js (r, t)= dv vf s (r, v, t) , 



(77) 



and then, knowing the concentration gradient, the diffu- 
sion coefficient from the macroscopic equation (jsj- 

Solving the linear Boltzmann-Lorentz equation (|75|l is 
by far not trivial, therefore, the Chapman-Enskog ap- 
proach has been developed to solve this equation approx- 
imatively. This approach is based on two simplifying as- 
sumptions: (i) f s (r, v, t) depends on space and time only 
trough the macroscopic fields and (ii) f s (f, v, t) can be 
expanded in terms of the field gradients. In the gradient 
expansion 

f s (r, v, t) = (f, v, t) + A/M (r, v, t) + A 2 /j 2) (?, v,t)... 

(78) 



a formal parameter A is introduced, which indicates the 
power of the field gradient. At the end of the computa- 
tion we set A = 1. 

The first term at the right-hand side, /i° (r, v, i), is 
the distribution function of the uniform system. Since 
the tagged particles are mechanically identical to the rest 
of the particles, fs°\r,v,t) is just proportional to the 
distribution function of the embedding gas: 



(79) 



Substituting Eq. i|78|) into the Boltzmann-Lorentz equa- 
tion l|75|l and collecting terms of the same order of A (that 
is of the same order in the gradients) we obtain successive 
equations for /j ' (r, v, t) , fs (?, v, t), etc. The zeroth- 
order equation in the gradients yields /i "* (r, v, i) for the 
homogeneous cooling state in accordance with Eq. I|79|l . 
The first-order equation reads 



d (o) /s (D 3C1) /s (o) 



dt 



dt 



■v-Vf 



(0) 



.92 



(a)/ (/,/«) , (80) 



where f /dt indicates that only terms of fc-th order 
with respect to the field gradients are taken into account 
when the time derivative of some function / is computed. 
For the case of interest the macroscopic fields are n s (f, t) 
and T (r , t) , which satisfy in the homogeneous cooling 
state the equations 



dn s fdW 9« 

~dt ~ \ dt + dt + ~dt 

dT_(d^ 0(i) 0(2) 

~dt _ \~~dT + ~&T + ~dt 



Comparing the terms in A we find 
0(o) n 9 (o) T 



n, = A 2 L>V 2 n, 



T = -cr. 



dt 



0. 



dt 



-CT, 



(81) 



(82) 



which yields for the first term in the left-hand side of Eq. 



0(°)/ s (1) d^n s dfP , 0(°)T dfi l) _ CT dfP 



dt dt dn s dt 

Similarly, taking into account 



dT 



dT 



d {1 



dt 



= 0. 



dt 



= 0. 



(83) 



(84) 



according to Eq. (|81|l . we obtain the second term in the 
left-hand side of Eq. (|SUj) . 



0W/i O) d^n s 0/ s (O) 0«T df ; 



(0) 



dt 



dt dn s 



+ 



dt dT 



(85) 
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Using Eq. I|79|l with constant n and /, we write the last 
term in the left-hand side of Eq. I|80l) , 



AO) 
p(0) _ J< 



v- V/ S w = ^v-Vn s = -tv-Vn s )f 



(86) 



and summing over i = j yields 



D= -- I dvv-G. 



(92) 



and obtain finally the equation for /J : 

C T ^+52(^(/,/i 1) ) =±(v-Vn s )f. (87) 

We search for the solution of Eq. H87fl in the form 

fP = G(v)-Vn s (?,t) , (88) 
which implies with Eq. I|77|) the diffusion flux 

J s = j v(G (v) ■ Vn s ) dv = -D\7n s , (89) 

where we take into account 

u/j°W=0 (90) 

due to the isotropy of fs°\\v |). From Eq. (|59^l follows 



J ViGjdv = —SijD 



(91) 



We substitute / s (1) from Eq. (|SS|l into Eq. ljS7|) discarding 
the factor Vn s . Then we multiply it by and integrate 
over v to obtain 



The first term in the left-hand side equals —(TdD/dT, 
while the right-hand side equals T/m, according to Eq. 
(|2*oTl . The structure of Eq. ljS7|) suggests the Ansatz 



G (v) oc vf(v) = bvf(v) . 



(94) 



where the constant b is to be determined from the above 
equation. With this Ansatz the second term in the left- 
hand side of Eq. I|93|) may be written as 



g2 g°^ / cfr/iffil (j, (f) = 92 ^ ^ / dvidv 2 f deQ (-v 12 ■ e)\v 12 ■ e\ f (v!) f (v 2 ) x (vx - v 2 ) ■ (vl - vx) 



= --^n 2 a 2 g 2 (<j)v T / dcidc 2 / de6 (-ci 2 • e) |ci 2 • e | / (ci) / (c 2 ) (1+e) (c 12 ■ e f = t, 

bm . m 



-l 

ad 



To evaluate the last expression we use the main property of the collision integral Q , 

dvx^J (tTi) /(/,/) = y / d«id« 2 / de (-ui2 • e ) |« 12 • e | / (p u t) f (v 2 , t) A [if, (ft) + $ (v 2 )\ , 



r 



(95) 



(96) 



the relation 



1+e 



v 2 - v 2 = - {v[ - vi) = (vi 2 ■ e ) e 



(97) 



which follows from the collision law Eq. Q and the sym- 
metry of the expression with respect to the exchange of 
the particles indices 1 *-> 2. In Eq. (|95|l we also use 
the scaled distribution functions / (ci) and /(c 2 ) and 
the definition of the adiabatic relaxation time T„ )a <i, Eq. 
(1711) . Taking finally into account that 



D — — — j dvv ■ vf(v) 



bTn 
m 



according to Eqs. lj9Hl9"I|) and we recast Eq. (jHSJ 



into the form 



(99) 



Simple arguments show that the diffusion coefficient 
scales as D cx VT: According to Eqs. f7T|l and 
the quantities r v ae j and £T scale as 



CT^; oc C ccvt oc VT. 



(100) 



At the same time, the right-hand side of Eq. (|99|l scales 
as oc T, which implies that D oc yT. Therefore 

dD D , . 

T 0T = J ^ 
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and we obtain the solution of Eq. (|99|l , 



T 

D(t) = - 
m 



(102) 



which coincides with Eq. 167|l . Substituting T„. a d and 
C given by Eqs. (JZU and (32) into Eq. (|102jl we ob- 
tain again the diffusion coefficient, Eq. (|73|1 . Hence the 
Chapman-Enskog approach yields for the case e =const. 
the same result for the diffusion coefficient as the Green- 
Kubo approach. 



VIII. DIFFUSION COEFFICIENT FOR 
IMPACT- VELOCITY DEPENDENT 
COEFFICIENT OF RESTITUTION 

A. Dynamics of a granular gas with 
impact-velocity dependent coefficient of restitution 

In the previous sections it was assumed that the coef- 
ficient of restitution, which determines the collision dy- 
namics in granular gases, is constant, i.e., independent of 
the impact velocity. Now we consider collision models, 
where this dependence is taken into account. 

For both models addressed above, the stepwise model 
due to Eq. © and the realistic dependence for viscoelas- 
tic spheres, Eq. @, there exist a characteristic velocity, 
namely g* for the former case and 7 -5 for the latter. The 
existence of this additional characteristic velocity, which 
is not related to the thermal velocity, leads to serious 
consequences for the granular gas dynamics: Rescaling 
time according to the Eq. (|55j) does not lead to an ef- 
fective steady-state of the system. Indeed, although the 
average kinetic energy of particles and their collision fre- 
quency is kept constant under this transformation, the 
value of the characteristic velocity (g* or 7 -5 ) in the 
new units increases with time. This follows from the fact 
that in the homogeneous cooling state the collision fre- 
quency steadily decreases in the laboratory time. Hence 
the time, measured in the new units, slows down [see Eq. 
(|53|l ] and the characteristic velocity (fixed in the labora- 
tory time) respectively grows. Therefore, a steady-state 
may not be achieved under this transformation, since the 
characteristic velocity, which determines the collision dy- 
namics, is not a constant in the new time units. 

Consequently, in contrast to the case e =const, it is 
not possible to map the dynamics of a granular gas with 
impact-velocity dependent coefficient of restitution onto 
a steady state dynamics of a system whose time slows 
down. As the result we cannot use Doob's theorem [2j| 
and deduce the exponential correlation function under 
the assumption of a Markovian collision process. This 
makes application of the Green-Kubo approach much 
more complicated, unless the adiabatic approximation is 
adopted The adiabatic approximation assumes that 
the relaxation rate of temperature, proportional to 
Eq. (32Jl, is negligibly small as compared to the relax- 
ation rate of the velocity distribution function, T v (t) . 



Under this approximation one can still assume exponen- 
tial form of the velocity correlation function, since on 
the time scale of the velocity correlation, the system is 
almost in a steady state |34j. 

From the above Eqs. (J3SJ, (33|l . lj4"5j) for £ and Eqs. 
(|62[1 . (|72|l for r„ follows for the case e = const.: £(i) ~ 
(1 — e 2 )r c (i) _1 and ^(t)' 1 ~ r c (i) -1 > which yields the 
condition of the validity of the adiabatic approximation, 



(1 - e 2 ) « 1 



(103) 



The same condition may be obtained for the case of 
impact-velocity dependent coefficient of restitution. If 
this condition is not justified it is preferable to use the 
Chapman-Enskog approach to compute the diffusion co- 
efficient. 



B. Diffusion coefficient for the stepwise 
impact-velocity dependent coefficient of restitution 

Since the shape of the velocity distribution function 
does not depend on time for this collision model [26| 
the derivation of the diffusion coefficient by means of 
the Chapman-Enskog method is identical to the case 
e =const, presented in the previous chapter. It leads 
thus, to the same Eq. I|99|l for the diffusion coefficient. 
The cooling coefficient is given by Eqs. (|35l I47f) for this 
model, while T~\ d may be found from the general rela- 
tion, Eq. jni- Using /(ci) = cj>(c{) and / (c 2 ) = <j){c 2 ) 
in Eq. I|71|) and Eq. @ for the coefficient of restitution, 
we again express T~^ d as a kinetic integral. Its solution 
yields (see Appendix lB|l 



'2 e 



- 1 



mg 



AT 



exp 



mg 



AT 



(104) 



with T- l (t) given by Eq. I|54[) . To find the diffusion 
coefficient, the first-order non-homogeneous differential 
equation l|99|) must be solved, which may be always writ- 
ten in quadratures. The resulting expression is rather 
complicated and not of practical use. In both limits 
T > mg* 2 /A and T <C mg* 2 /A the diffusion coefficient 
scales in the same way as for e =const, that is, D oc \/T. 
Therefore in these limits Eq. I|102[) stays valid. We ex- 
trapolate this dependence for the full interval of temper- 
ature, that is, we use the approximation 



D{t) 



T 



1 



v, ad 



(<) - oCW 



(105) 



with £ and t v , a d given respectively by Eqs. (|35l I47J1 and 

The interesting feature of this model is the relatively 
sharp dependence of the diffusion coefficient on the gran- 
ular temperature: For T mg* 2 /A the diffusion coef- 
ficients varies exponentially fast with temperature. If 
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e* = 0.6 it changes by about 50% in a narrow interval 
of temperature, Fig. ^ Note that the approximation 
(I105J1 agrees fairly well with the result obtained by the 
numerical solution of Eq. (|99|l . Fig. 2] 
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FIG. 1: Reduced self-diffusion coefficient D/Dq, where Do is 
the Enskog self-diffusion coefficient given by Eq. 1741 . as a 
function of temperature for the stepwise impact-velocity de- 
pendent coefficient of restitution, Eq. @. Full line: numeri- 
cal solution of Eq. 1)990 : dashed line: approximation given by 
Eq. 11021 . The parameters are: m — 1, a = 1, g* = 1. 



C. Diffusion coefficient for gases of viscoelastic 
particles 



The application of the Chapman-Enskog method for 
granular gases of viscoelastic particles is similar to the 
case e =const. For gases of viscoelastic particles, how- 
ever, some complications arise which are related to the 
additional dependence of the velocity distribution func- 
tion on time: Here the shape of the velocity distribution 
function is not invariant, but depends on time via the 
time-dependent coefficients of the Sonine polynomial ex- 
pansion 0, l4o| . If we keep only two first non- vanishing 
terms of the Sonine polynomial expansion for the distri- 
bution function, we need to account for the time depen- 
dence of the second Sonine coefficient a 2 (t). Then Eq. 
(|83|l adopts the more general form 

9(0) fW d (o) ns df W g(0) T g f W g(o) a2 df W 



dt 



dt dn s 



dt dT 

dfP 
dT 



a 2 ■ 



dt 
da 2 



da 2 



(106) 



If we again choose /, 



as given by Eq. (|SS|l and 



(j5fj|l to Eq. (|!2III), we arrive at 



1 

3° 2 



dv', 



dG 
da 2 



92 (cr) 



3nm 



mv 2 
dv j [V 



(107) 



We substitute G(v ) by means of Eq. (|MJ> with 

/(«, t) = f M {v) [1 + a 2 {t)S 2 (c 2 )] , (108) 

where Jm(v) = (n/u|>)0(c) is the unsealed Maxwell dis- 
tribution. 

All terms in Eq. (|107|l . except for the second term 
in the left hand side, have been already evaluated. The 
remaining term reads 



a 2 f dG a 2 

— / dv v ■ = — b 

3 J da 2 3 

= -a 2 b nv T 



dv v ■ vf M S 2 (c 2 ) 

)c 2 S 2 (c 2 )=0, (109) 



dc 



where we use 



dc(j)(c)c 2 S 2 (c 2 ) = , 



(110) 



which follows directly from the definition of the second 
Sonine polynomial, Eq. 139fl . Hence we arrive again at 
Eq. l(9"9"|) for the diffusion coefficient as for the case e = 
const: 



(111) 



Evaluating the kinetic i nteg ral, Eq. (|71|l . for viscoelastic 
particles yields t"^ [Hilf: 

2 



-1 

. ad 



3 Tc 



l (t) 



1 + r^-7%= 6 '^ 

16 4V27T 



(112) 



where r c 1 (t) is given by Eq. (|54|l . while a 2 {t) and 8'(t) 
by Eqs. 1491 150|) . Hence we note, that both T~* d as well 
as C with Eqs. i|35l 15 1|> may be written as an expansion 
in terms of the small parameter 5' . 

To solve Eq. I|lll|) for the diffusion coefficient we write 
D(t) as an expansion too, 

D{t) =D (l + 5'b 1 +5' 2 b 2 + ---^ (113) 

and substitute this expression into Eq. (|lllf> together 
with equivalent expansions for ((t) and T~\ d (t). The 
resulting equation can be solved perturbatively for each 
order of 6'. Here we give the final result of these straight- 
forward calculations 0, Eal : 



Pit) _ \ 1 169 co 
Do [ 320 

( 2867777 lu 2 3611 w x 



perform the same steps of analysis which lead from Eq. 



V 20480000 7T 7040 ^ 



8 I2 + 



(114) 
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where Dq is the Enskog coefficient of self-diffusion at the 
initial temperature T , 



O = r52 (a)W — 



(115) 



IX. CONCLUSION 

We studied diffusion of tracer particles (self-diffusion) 
in a gas of dissipatively colliding particles. Two main 
theoretical approaches for the calculation of transport 
coefficients have been considered: the Green-Kubo and 
Chapman-Enskog approaches. The Green-Kubo ap- 
proach is based on the computation of the time integral 
of the time-correlation functions of dynamical variables, 
while Chapman-Enskog method is based on the Boltz- 
mann equation. The dissipative collisions of the parti- 
cles are described by the coefficient of restitution. We 
addressed three different models for this coefficient: (i) 
the model of a constant coefficient of restitution, (ii) the 
model of a simplified (stepwise) dependence of the coef- 
ficient of restitution on the impact velocity, and (iii) the 
realistic model for the coefficient of restitution which fol- 
lows from the contact mechanics of viscoelastic particles. 

For the case (i), e =const, the dynamics of the dissi- 
pative system may be mapped onto the steady-state dy- 
namics of another gas. This allows for the application of 
the standard methods of statistical physics, elaborated 
for molecular gases in thermodynamic equilibrium. In 
particular the velocity time-correlation function may be 
straightforwardly evaluated and the diffusion coefficient 
is correspondingly found. The obtained Green-Kubo dif- 
fusion coefficient coincides with that derived using the 
Chapman-Enskog approach. In contrast, for the case 
of impact-velocity dependent coefficients of restitution, 
models (ii) and (iii), there is no simple transformation 
to map the dynamics of the system onto an effective 
steady-state dynamics. This turns the application of the 
Green-Kubo method either restrictive (e.g. by assuming 
an adiabatic steady-state), or very complicated and in- 
efficient. Therefore, for gases of particles which interact 
by an impact-velocity dependent coefficient of restitu- 
tion, we use the Chapman-Enskog approach to compute 
the diffusion coefficient. 

Although the Green-Kubo approach is more general, 
in principle, since it does not assume any form of the 
time correlation function, while the Chapman-Enskog ap- 
proach exploits the assumption of molecular chaos, the 
Chapman-Enskog method appears to be preferable for 
the application to granular gases. While the latter one is 
applicable to any collision model, the former one becomes 
either restrictive or very complicated if the coefficient of 
restitution depends on the impact velocity. 



APPENDIX A: SECOND SONINE COEFFICIENT 

To find the velocity distribution function in the homo- 
geneous cooling state we write f(y) in the rescaled form, 
Eq. (|33|l . and substitute it into the Boltzmann equation 
(|75|l (with the term v ■ V omitted and with f s substituted 
by /). Then this equation may be reduced to two equa- 
tions, Eq. (|34|l for temperature and another one for the 
scaled distribution function / (c, t) 0, : 



|( 3 + Cl ^)/ et) + B-.|/(c- t ) 



I 



where 



B = B(t) = v T (t)g 2 (a)a 2 n 



(Al) 



(A2) 



and 7 (f, is the dimensionless collision integral, that 

is, the collision integral Eq. I|75|l written in terms of the 
scaled distribution functions and scaled velocities. The 
analysis shows that on the time scale of the diffusion pro- 
cess - the hydrodynamic time scale - the term _B -1 ^-/ 
in Eq. IjAlf) can be omitted |45| . 

The solution for the scaled distribution function / may 
be found as an expansion around the Maxwell distribu- 
tion in terms of orthogonal Sonine polynomials S p (x) 



m = 0( C ) 



l + 5> p (i)S P (c 2 ) 



(A3) 



where Sq(x) — 1, S\(x) = 3/2 — x and S 2 (x) is given by 
Eq. gH. 

Multiplying both sides of Eq. IjAlf) by cf , integrating 
over c\, and using the orthogonality relation for the So- 
nine polynomials, we obtain an infinite set of equations 
for the moments fip 0, 03 



3 Hp = (J.2 V (c p ) 



p = 2,4,- 



(A4) 



Here (c p ) = J dcc p f(c) are the moments of the velocity 
distribution function, 



jip = — J c p Idc 

= ~ J dcidc 2 J deQ (-cj. 2 ■ e) |c*i2 ■ e\ 

x f( Cl )f(c 2 )A(c p + c p 2 ), (A5) 

and we use the main property of the collision integral, 
Eq. H96JI . Since (c p ) and fj, p are expressed in terms of the 
Sonine coefficients, the set of equations ljA4() may be used 
to determine these coefficients, i.e., to find the velocity 
distribution function. Assuming a cutoff of the series 
Eq. (|A3|) . i.e. assuming that the Sonine coefficients a/, 
with k > kg are negligible, the infinite set of equations 



14 



(|A4|) turns into a closed finite set of equations for the 
coefficients a^. 

The first Sonine coefficient oi = 0, according to the 
definition of temperature [i(J, Since (c 2 ) = 3/2, 

the first equation in Eq. I|A4(1 for p = 2 is an identity. 
Assume that 03, 04, . . . may be neglected. Then using 
(c 4 ) = (15/4) (1 + 02) which follows from the definition 
of the second Sonine polynomial, Eq. I|41|l . and the ap- 
proximation /(c, t) — 4>{c)[l + a2(t)S l 2(c 2 )], we write Eq. 
(|X4|) for p = 4: 

5/i 2 (l + a 2 )-/i 4 = 0. (A6) 

The coefficients /x 2 and p^ are functions of 02 and have 
the form of kinetic integrals, Eq. (|42(l and may be 
straightforwardly evaluated by the scheme introduced in 
For the case e — const., the coefficient /X2 is given 
by Eq. gSJ), while // 4 reads [HUEI 

Pi = 4^MTi + a 2 T 2 ) (A7) 

with 

4 V 2 / ( A8 ) 

T2= ll (l ~ £2) (69 + 10e 2 ) + i(l + £ ). 

Substituting Eqs. I)A7|) and (g^) into Eq. (|ABf> we ob- 
tain an equation for the second Sonine coefficient a 2 . Its 
linear solution is given by Eq. (1441) . while the next-order 
solution is presented in Refs. |1. l42l|. 

The velocity distribution function for a gas of viscoelas- 
tic particles may be found just in the same way: We eval- 
uate H2 and ^ which have the form of kinetic integrals, 
substitute them into Eq. (|A6J) and solve the resulting 



equation for the second Sonine coefficient 0,2- The final 
result is given in Eq. 149|) . Details of the derivation can 
be found in |3.l4fij. 

The derivation of the Sonine coefficients, outlined 
above, is based on the hypothesis of molecular chaos 
which is the main precondition of the Boltzmann equa- 
tion. In the coarse of the evolution of the granular gas, 
correlations between the particles which have been ne- 
glected in this approach, may become important. In par- 
ticular, it was found in MD simulations |44j that accu- 
mulating correlations lead to deviations of 02, from the 
theoretical predictions, e.g. [ioj . 



APPENDIX B: KINETIC INTEGRALS FOR 
STEPWISE IMPACT- VELOCITY DEPENDENT 
COEFFICIENT OF RESTITUTION 

The scheme elaborated in 0,^1 f° r the computation of 
kinetic integrals may not be directly applied for the eval- 
uation of the cooling coefficient £ and the relaxation time 
t v , ad in case of the stepwise impact- velocity dependent 
coefficient of restitution defined in Eq. (0. Therefore, 
here we outline their evaluation. 

By means of Molecular Dynamics simulations it has 
been shown that in a granular gas with coefficient of resti- 
tution given by Eq. the second Sonine coefficient is 
negligible, that is, 02 ~ 0, provided e* < 1. [26^. Then 
with 

(1 + e) = 2 + (e* - 1)6 [-g* - (c 12 • e )] , (Bl) 

where g* = g* /vt, we write the factor in t~ j(t), which 
has the form of the kinetic integral as 



2 J dcidc 2 J deQ (-c 12 ■ e) \c X2 ■ e \ f (ci ) / (c 2 ) (ci 2 • e f 

+ {e* -I) J dc x dc 2 J deQ (-c 12 • e ) |c 12 • e | / (ci ) / (c 2 ) (c 12 • e f 6 [-g* - (c 12 • e )] . (B2) 



Changing variables, (ci, c 2 ) -> [C, c 12 j , where C = (c x + Substituting Eq. $52$ with Eq. $5$ into Eq. JB2J) and 
c 2 )/2 and C12 = ci - c 2 , we can write takin S into acc °unt that J <p{C)dC = 1, we obtain for 



/(cx)/(c 2 ) = 0( Cl2 )^(C) (B3) 

with 



the second term in Eq. (|B2|1 

/>oo 

(e*-l)47r/ c\ 2 (j)(ci2)dci 

J a* 



( Cl2 ) = (27r)- d/2 exp(-ic 2 2 l x / sin 6 ^ cos " 6 d0 I ( B5 ) 



( , V I2 (B4) 
<KC) = f - J cxp (-2C 2 ) 



2 iz / "'o Jo 



2- 



2^ ^ ^ r > r<2\ with O = g*/ci 2 . In Eq. (|B5|) we take into account 

that the integral vanishes if the normal component of 



15 



the impact velocity \c\2 • e\ — |ci2COS#| is smaller than 
g*. Simple calculations then yield for the second term in 
Eq. (jB2j: 

^<e--l)(l + ^)«p(-^). (B6) 

The first term of the kinetic integral in Eq. I|B2|) may be 
obtained from the second term, taking g* — and omit- 
ting the factor (e* — 1). Summing up these two terms 
and multiplying the result by (1 / '6)vt92{<^)<^ 2 n, accord- 
ing to the definition Eq. (|71|) . we arrive at Eq. (|l()4f> for 



The cooling coefficient £ may be evaluated in a similar 
way. From the definition of the collision model, Eq. 
and the general collision law, Eq. Q), follows 

A (c\ + c|) = - (1 - e 2 ) (c 12 • e f 

= - (l - e* 2 ) (c 12 • ef 6 [-g* - (c 12 • e )] . (B7) 

Substituting the latter expression into Eq. HlffiJ) for /i2j we 
arrive at the same kinetic integral as for T~* d (i), which 
eventually leads to the final result Eq. Ij47(l . 
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